Mon. Not. R. Astron. Soc. 000,[T]{8](2012) Printed 9 October 2012 (MN ET^X style file v2.2) 



An empirical prediction for stellar metallicity distributions in 
nearby galaxies 

Molly S. Peeples 1 * and Rachel S. Somerville 2 

1 Southern California Center for Galaxy Evolution Fellow, University of California Los Angeles, 430 Portola Plaza, Los Angeles, CA 90095 

2 Department of Physics and Astronomy, Rutgers, The State University of New Jersey, 136 Frelinghuysen Road, Piscataway, NJ 08854 



(N 

o 

-i— > 
O 

O 

oo 



o 

U 

6 



(N 
> 

o 



l October 2012 



ABSTRACT 

We combine star-formation histories derived from observations of high redshift galaxies with 
measurements of the z ~ relation between gas-phase metallicity, stellar mass, and star 
formation rate to make an explicit and completely empirical connection between near-field 
and distant galaxy observations. Our approach relies on two basic assumptions: 1) galaxies' 
average paths through time in stellar mass vs. star formation rate space are represented by 
a family of smooth functions that are determined by the galaxies' final stellar mass, and 2) 
galaxies grow and become enriched with heavy elements such that they always evolve along 
the mass-metallicity-star formation rate relation. By integrating over these paths, we can 
track the chemical evolution of stars in galaxies in a model independent way, without the 
need for explicit assumptions about gas inflow, outflow, or star formation efficiency. Using 
this approach, we present predictions of stellar metallicity (i.e., O/H) distribution functions 
for present day star-forming galaxies of different stellar masses and the evolution of the a- 
element stellar metallicity-mass relation since z <~ 1. The metallicity distribution functions 
are fairly well described as Gaussians, truncated at high metallicity, with power-law tails to 
low metallicity. We find that the stellar metallicity distribution for Milky Way mass galaxies is 
in reasonable agreement with observations for our Galaxy, and that the predicted stellar mass 
vs. mean stellar metallicity relation at z = agrees quite well with results derived from galaxy 
surveys. This validates the assumptions that are implicit in our simple approach. Upcoming 
observations will further test these assumptions and their range of validity, by measuring the 
mean stellar mass-metallicity relation up to z ~ 1, and by measuring the stellar metallicity 
distributions over a range of galaxy masses. 

Key words: stars: abundances — ISM: abundances — Galaxy: abundances — Galaxy: stellar 
content — galaxies: abundances — galaxies: evolution — galaxies: stellar content 
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1 INTRODUCTION 

As stars are born with the imprint of the elemental abundances in 
their gaseous birth environments, measurements of stellar metal- 
licities within galaxies can be us ed to gain ins i ght into their star 
formation and chemical histories i Tinslev|[l975l : Tinslev & Larson! 
ll978l : lBland-Hawthorn etaitololT In this paper, we carry out an 
explicit test of the self-consistency between empirically derived 
star formation histories, observational scaling relations between gas 
phase metallicity (Z g ), stellar mass (M*), and star formation rate 
(Mspr), and observed stellar metallicities at z = 0. 

It has long been known that there is a strong positive correla- 
tion between the luminosity or stellar mass of st ar forming galaxies 
and their galaxy-averaged gas-phase abundance dGarnett & Shields] 
ll987l : IZaritskv et al.ll994l) . This relationship (the "mass-metallicity 
relation") has now been measured for large samples of nearby 
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galaxies (e.g. iTremonti et al.l |2004 iKobulnickv & Kewlevl |2004 
iKewlev & Ellison! 2008)) and for smaller samples of distant galax- 
ies out to z ~ 3^5 llSavaglio et al.l 120051; IShaplev et all 120051 : 
lErb et alj|200rj : iMaiolino et alj|2008l ; IMannucci et alj|2009h . The 
normalization, slope, scatter, and evolution of the observed mass- 
metallicity relation place strong constraints on theoretical mod- 
els of galaxy formation, in particular on the physics of feedback 
from massive stars and supernovae. Recently, it has been sug- 
gested that the scatter in the observed mass-metallicity relation 
can be red uced by considering the star formation rate a s a third 
parameter dLara-Lopez et alj|201fj : IMannucci et al. Physi- 
cally, the star formation rate is probably acting as a proxy for the 
galaxy gas fractions, as galaxies with higher gas fracti ons have both 
highe r star format ion rates and more dil uted metals dHug hes et al .1 
l20ll . Moreover, lLara-Lopez et al] and IMannucci et al] find that 
galax ies up to z ~ 2.5 (e.g., ISavaglio et al.l 120051 : lErb et al .1 
1 20061) appear to always lie on the same M t -Z g -MsFB. "funda- 
mental" relation, though this may begin to break down at z ~ 3 
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dMaiolino et~ai]|2008l . though see also lYabe et~al]|2012h . It is not 
yet clear whether or why a non-evolving mass-metallicity-star- 
formation rate relation should be genericall y expected , altho ugh 
some cosmological models do predict this foavg et alj|2012l M. 
Arrigoni et al. in preparation). Self-regulati on due to stellar-driven 
outflows is u s ually i nvoked in this context jFinlator & DaveJ | 2008 |; 
Spitoni et alj 1201(1 IPeeples & Shankari 1201 ll ; iDave et alj hoill ; 



Daval et alj|2012l) 



There is a long history in the literature of attempts to model 
the chemical evolution of galaxies using a "classical" approach, 
the simplest version of which is the "closed box", in which galax- 
ies start out with all of their mass in the form of pristine gas, 
and convert some fraction of that gas to stars, producing a certain 
"yield" of heavy elements along the way. It was quickly realized 
that the simplest closed box picture could not reproduce observa- 
tions, in particular of the distribution of stellar metallicities in the 
Solar Neighborhood in our Galaxy (the "G-dwarf" problem), lead- 
ing to modifications such as inflowing p ristine or pre-enriched gas 
as well as outflows of ga s and metals dLvnden-Belllll975l ; IPagell 
ll989tlColavittietalj2008l) . While much has been learned from this 
approach, it has the obvious drawback that the results tend to be de- 
pendent on the rather arbitrary functions adopted for the inflow and 
outflow rates. 

Another approach that has become fairly widely used is 
the coupling of detailed chemical evolution models with cos- 
mological models of galaxy formation, realized in the form of 



either numerical hydrodyna mic simulatio ns (|S cannapie co et al 



2006; Kobavashi et al. 2007; Brooks et al. 2007; Mouhcine et al 



2008; Oppenheimer & Dave 2008; Wiersmaet al. 2009) or semi 



analytic models l De Lucia et alj 20041; iNagashima et alj 1 20051 ; 
ISomerville et alj2008l ; lArrigoni et alj20ld)~ These have the advan- 
tage that quantities such as the rate of inflow of gas into galaxies is 
motivated by the formation rate of structure in the Cold Dark Matter 
(CDM) paradigm that provides the backbone for these models, and 
outflow rates are also governed by physically motivated scalings. 
However, they suffer from the drawback that currently, all CDM- 
based galaxy formation simulations, whether numerical or semi- 
analytic, apparently fail to reproduce the observed scaling of star 
formation history functional for m with galaxy mass ( sometimes re- 
ferred to as "downsizing"; see iFontanot et atll2009l for a detailed 
discussion). 

Here we propose a different approach to modeling the buildup 
of metals in galaxies that sidesteps some of the problems with 
both the classical and cosmological methods. We make use of a 
family of parameterized empirical star formation histories, derived 
from observations of th e stellar mass vs . star formation rate relation 
at different redshifts bv lLeitnerl ( f2012h . These empirically-derived 
star formation histories are in good agreement with those derived 
from spectral energy distri bution modeling and the evolution of the 
stella r -halo mass relation llZheng et alj|2007 | ; [ Conrov & WechsleJ 
120091 ; iBehroozi et alj|2012l ; iMoster et alj|2012l) . We assume that 
the average star formation history for galaxies with a given stel- 
lar mass today can be represented by this smooth functional 
form. This assumption is supported by mounting observational ev- 
idence that galaxies with masses of roughly the Milky Way or 
lower (~ few xlO lo M0) build up most of their mass through 
smooth accretion of g as, with mergers and starbursts playing 
a relatively minor role jNoeske et alj|2007l; iRobaina et alj 120091 ; 
lOliver et alj|2010l ; iKarim et alj|201 ll ; lLeitnerll2012l) . This is also 
supported by results from cosmological simulati ons l lBrooks et alj 
l2009l ; lHirschmann et alj2"012l ; lTissera et al.l2012h . Then, once stel- 
lar mass loss and recycling has been accurately accounted for 



dLeitner & Kravtsovll201 ll) . the growth of galaxies can be traced 
through cosmic time by stepping along these Msfr(M*) relations 
as a function of redshift. 

We make a second assumption, that galaxies always evolve 
along the M^-Zg-MsFn relation. That means that as galaxies 
evolve along their trajectories in M*-Msfr space (specified by 
their star formation history), we assume that the stars born at that 
moment form out of gas with a metallicity Z g given by the M*- 
Z g -MsFR relation. In this way, we can build up the distribution 
function of stellar metallicities for galaxies with different masses 
at z = 0. We can also test whether the two assumptions of em- 
pirical star formation histories plus a non-evolving M*-Z g -MsFR 
relation lead to an average stellar metallicity relation at z — that 
agrees with observations. Both of these quantities will be better 
constrained observationally over a wider range of galaxy masses, 
and to higher redshift, in the near future. 

This paper is organized as follows. In §|2] we describe our 
adopted and predicted scaling relations connecting galaxy stellar 
masses, star formation rates, and gas-phase metallicities across cos- 
mic time. In §[3]we present the stellar metallicity distribution func- 
tions and evolution of the stellar metallicity-mass relation predicted 
by this model. We conclude with a brief discussion of relevant im- 
plications in §|4] 

Throughout we assume a IChabrieJ J2003h stellar initial mass 
function (integrated o ver 0.1-100 M w ) and 12+log(0/H) = 8.7 
jAsplund etal1l2009l) . although we note that the 8.9 value found 
from helioseismology ( Delahave & Pinsonneauli 2006T) wo uld lead 
to 0.2 dex lower Z/Zq values than stated here. Following iLeitne j 
< |2012|) . we adopt a flat ACDM cosmology with Qm = 1 — = 
0.258 and h = 0.72. 



2 EVOLUTION OF SCALING RELATIONS 

The model we present here is based on two simple assumptions: 1) 
galaxies' average star for mation histori es are given by the empiri- 
cal functions provided bv lLeitneil d2012h . in which their z = stel- 
lar mass determines the shape of their star formation histories. 2) 
Galaxies evolve along a "fundamental" M*-Z g -MsFn relation that 
is itself unchanging with redshift. Taken together, these pieces pro- 
vide a unique enrichment history for an average galaxy of a given 
present day stellar mass M*,o. In this section we describe these two 
pieces of our model, where they come from, and some caveats, in a 
bit more detail. 

The starting point for the first piece of the model is mea- 
surements of g alaxy stellar masses and star formation rates from 
z ~ to 2.5 bv lKarim et all d201lh . from the COSMOS survey. It 
is now well established that M* and A/sfr show a fairly tight and 
nearly linear relationship which evolves in nor malization, but very 
little in slope, over a very broad redshift range jNoeske et al.l2007l ; 
iDaddi et alj|2007l ; IKarim etal.ll201ll ; lElbaz et alj|201ll) . Thesere- 
sults of course represent snapshots of d ifferent galaxy populations 
at different cosmic times. lLeitnejj2012l) first fit power laws to these 
snapshots in time, using the Msfr(M+, 2) results measured by 
iKarimetalJ l feOllh , 



M SF r(M*,z) = An 



10 11 M f 



0+1 



(1 



(1) 



where An = 0.0324 x 10 11 M Q Gyr~ 1 , a = 3.45, and /3 = 
—0.35. He then characterized how individual galaxies would have 
to grow over time in M* vs. Msfr such that the evolution of these 
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Figure 1. Evolution of the Af*-AfsFR relation as a function of redshift. Colored lines show "snapshots" of the galaxy population at different redshifts; growth 
of individual galaxies is shown by black curves with dots. Top right: Evolution of the mass-metallicity relation as predicted by the M*-Z g -AfgpR relation. 
Bottom left: Evolution of the Z g -AfsFR relation. Bottom right: Evolution of Z s as a function of the specific star formation rate, Msfr/M*. In all panels, 
individual galaxies are shown in black, from left to right, with log AI^/Mq = 8, 9, 10, 11 at z = 0. 



relations would be satisfied, accounting for mass loss and recy- 
cl ing from evolving stellar po pulations using the method presented 
in lLeitner & Kravtsovl d20 1 lh . Thus, these Msfr(M*, z) relations 
can be integrated to give the star-formation history for a galaxy of 
a given stellar mass M,,o at z = 0, MsFR(iW*,n, £0 The top left 
panel of Figure[T] shows the fits representing the population aver- 
ages at different redshifts, with the trajectories of individual galax- 
ies of different Ai*,o overplotted. 

The two obvious limitations of this method are that it only 
describes galaxies that are actively forming stars at z = (i.e., 
there is no accounting for a galaxies that cease forming stars and 
"fall off" of these relations) and that it is only valid up to redshifts 
and stellar masses where the M*-Msfr relation is reliably mea- 
sured. However, the majority of galaxies with stellar masses in the 
range 10 9 < < 3 x 10 10 Mq are on the "star forming main 

1 A similar approach, with qualitatively similar results, has been presented 
bv lNoeske et alj2007l 



seque nce" terinchmann et alj2004l ; ISalim et alj 2007; iKimm et al.1 
120091) . i.e. on the relations used to derive the lLeitnea ( 120 1 2f) star for- 
mation histories. Bas ed on observational estimates of merger rates 
dRobaina et alj|2009h. th e scatter in the observed M ^-Msfr rela- 
tion jNoeske et all2007r) . and theoretical predictions dBrooks et al.1 
|200C galaxies in this mass range are expected to build up the ma- 
jority of their mass through smooth accretion, with mergers and 
starbursts playin g a mino r role. On th e second poin t, the star for- 
mation histories iLeitneJ derives from iKarim et alt s data are not 
entirely reliable at z > 1.6 due to observational uncertainties 
(i.e., the redshifts for which A/sfr has not been well-measured 
for M* < lO 9 ' 6 M Q ). We further note that the star formation his- 
tories for M* — 10 s M© galaxies at z = are entirely based on 
extrapolations, as the progenitors of these galaxies cannot currently 
be observed at higher redshift. This is true to some extent at some- 
what higher masses as well (10 8 < Af* < 3 x 10 9 M q), with 
extrapolation becoming less important as Af^o increases. iLeitnej 
notes that the star formation histories derived for A'h ~ 10 8 Mq 
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galaxies predict a much larger fraction of recent star format ion than 
those derived from resolved stellar population studies: the iLeitner] 
star formation histories for these low mass galaxies imply that all 
of their star formation takes place at z < 1, while the resolved stel- 
lar population studies indicate that a significant f raction of past star 
formation occurred at z 2> 2 l lWeisz et al .11201 lh . 

The other three panels show the predicted evolution of the Z g - 
M*, Zg-MsFR, and Z g -M*/MsFR relations, where for the gas- 
phase metallicity we plot the observable 12 + log(0/H) as would 
be measured from star-forming H II regions derived from the z ~ 
Af*-Zg-M S FR relation. While several fits to the M*-Z g -M S FR 
relation exist in the literature, we choose to calculate metallicities 
from the double-quadratic fit 

[12 + log(0/H)] FM R = 

8.90 + 0.37m - 0.14s - 0.19m 2 + 0.12ms - 0.054s 2 , ^ 

where m = log M* — 10 and s = log Msfr, fit t o Sloa n 
Digital Sky Su rvey (SDSS) data by IMannucci etHl fcoicb . 
iMannucci et al .1 postulate that this so-called "fundamental metal- 
licity relation" (FMR) holds at all redshifts, and show that the 
observed z < 2.5 galaxies fall along this relation when extrap- 
olated to higher star formation rates. Equation l[2j is fit over the 
ranges 9.1 < logM* < H-35 and -1.45 < M S fr < 0.80 (see 
IMannucci et alj|2010l for the exact Msfr [A/*] bins used in their 
calculation); we must extrapolate outside of this range in order to 
describe the metallicity evolution of 8 ^ log M A /Mp) ^ 11. 5 
galaxies at z = 0. This parameterization is not Man nucci et alj 's 
best fit to the data, but we find that this functional form is bet- 
ter suited to the necessary extrapolation than alternatives. Specifi- 
cally, Equation l(2j will under-predict the metallicities for low-mass 
galaxies, compared to a more robust fit using /10.32 = logM* — 
0.32 log Msfr as the independent variable, 

[12 + log(0/H)] FM R, M = 

(3) 

8.90 + 0.39s - 0.20s 2 - 0.077s 3 + 0.064s 4 , 

where s = /io.32 — 10 ( IMannucci et alj[2010l) . Unfortunately, this 
quartic parameterization leads to an up-turn in 12+log(0/H) when 
extrapolated to low /Lto.32; while no extrapolation in /io.32 from 2 = 
to z ~ 2 is needed to describe observed galaxies, extrapolation is 
required to describe the evolution of the loca l galaxies we wish to 
characterize. Finally, [Lara-Lopez et aljfeoiol) describe the relation 
between stellar mass, gas-phase metallicity, and star formation rate 
as a plane, 

[12 + log(O/H)] U0 = 

(4) 

0.891 log M* - 0.422 log M SF r + 0.086. 

We choose, however, to not adopt this formulation because it gives 
even lower metallicities at high-redshift than Equation (O, while 
overestimating metallicities at high M* . 

Finally, we note that uncertainties in the calibration meth- 
ods for determining 12 + log(0/H) compound the difficulty 
in accurately character izing the M*-Z g -AfsFR relation (e.g., 
iKewlev & Ellison! 120081) ; these uncertainties not only affect the 
slope of the relation, but also its normalization at the ~ 0.3 dex 
levelQ Thus it will certainly be possible to improve on the results 

2 Comparisons of gas-phase and stellar metallicities could provide an 
potentially useful avenue for placing combined constraints on "reason- 
able" ranges for the Solar oxygen abundance, normalization of the mass- 
metallicity relation, and the nucleosynethetic oxygen yield. 



presented here once the M*-Zg-MsFR relation has been measured 
to lower masses (e.g., B. Andrews & P. Martin i, in preparation) and 
with more robust metallicit y calibrations (e.g.. lNicholls et alj2012t 
lLopez-Sanchez et alj|2012l ; M. Dopita et al. in preparation). 

From the top-right panel of FigureQ] it is obvious that on short 
timescale s, galaxies evolve along the mean mass-m etallicity rela- 
tion (c.f. lPeeples & ShankaJboi ll ; |Pave et alj|2012l) . Though it is 
subtle, the mass-metallicity relation evolves slightly more rapidly 
at low masses than at high masses. The bottom two panels show 
that this is because massive galaxies asymptotically plateau in their 
metallicities, while low-mass galaxies are still rapidly evolving at 
late times. The galaxy tracks in the bottom-left panel of FigureQ] 
show the metallicities at which galaxies are at the peak of their star 
formation. In §[3] we integrate along these tracks to derive stellar 
metallicities as a function of cosmic time. 



3 AN EMPIRICAL PREDICTION FOR STELLAR 
METALLICITIES 

Because stars are born with the chemical imprint of the gas they 
formed out of, by tracing galaxies through the evolutionary tracks 
plotted in FigureQ] we are able to track — in an average sense — 
the buildup and change in stellar metallicities in different galax- 
ies. Figure|2] shows the distribution of the stellar [a/H] within in- 
dividual galaxies of different A/* at z — and the predicted rela- 
tion between galaxy-averaged metallicities (Z*) and Ah. Though 
the distribution o f stellar [Fe/H] ha s been well-measured for the 
Milky Way feg.. Ilvezic et al.|[2008h . [a/H] has not been as well 
characterized^ The histograms in the top panel of Figure|2] are 
based on all star formation, and thus are only applicable for 
comparing to metallicity distributions of low-mass stars, such as 
for the upcoming SDSS-III APO Galactic E volution Experiment 
(APOGEE) observations of the Milky Way dAllende Prieto et al] 
120081 : lEisensteinetai]|201ll) and M* ~ 10 s M Q dwarf Irregular 
gala xies in the Local Group (E. Kirby et al. in preparation; see 
also iKirbv et alj|20"l(i 1201 l j) . We note that our assumption that 
all galaxies of a given present-day stellar mass follow the same 
smooth star formation history may tend to underestimate the width 
of the ensemble-averaged stellar metallicity distribution function. 
It is well known that galaxies have somewhat stochastic star for- 
mation histories (although apparently they typically do not deviate 
from the star formi ng main sequence by more than a factor of 2-3; 
iNoeske et afll2007l) . and that the overall shape of the star forma- 
tion history can show considerable diver sity from galaxy to galaxy, 
particularly for low mass objects (e.g., IWeisz et al]|201ll) . These 
effects would tend to broaden the predicted metallicity distribution 
function and to lead to a corresponding galaxy-to-galaxy diversity 
in this quantity. However, they should not significantly affect our 
predictions for the galaxy- and ensemble-averaged mass metallic- 
ity scaling relations that we discuss below. 

The logM^/M© = 10.5 (bright g reen) distribution i s in- 
cluded for comparison to the Milky Way dKfypin et al.1l2002h : wc 
also sho w the distribution of [O/H] as measured for 72 local disk 
stars by l lBensbv et al]|2004l) . though we caution that it is unclear 
the extent to which these stars are representative of the entire Milky 
Way. The dashed portions denote where Equation (O has been ex- 
trapolated. At higher stellar masses, the metallicity distribution 



3 In an effort to keep our model as empirically-based as possible, we do 
not try to model Iron production directly. 
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log(a*/a Q ) AM/A log(OZH) 



-3.3907 
-2.9125 
-2.5897 
-2.3483 
-2.1570 
-1.9997 



0.3202 
0.3203 
0.3205 
0.3206 
0.3208 



-0.00006 
-0.00011 
-0.00024 
-0.00042 
-0.00066 
-0.00097 



2.23054 
2.22968 
2.22867 
2.22769 
2.22358 



Table 1. Stellar metallicity distribution function for M*,o = 10 10 5 M 
galaxy as plotted in the top panel of Figureff] Negative AM/ A log(0/H) 
denotes where Equation (2) has been extrapolated. Full tables for 
logM* i0 /M Q = 8, 9, 10, 10.5, 11 and 11.5 are available in the elec- 
tronic version. 



functions turn over; this is directly related to the decline in star 
formation rates at later times, resulting in the formation of rela- 
tively few higher-metallicity stars (see the bottom-left panel of Fig- 
ure!]}. Because the turnover in the star formation histories is rela- 
tively slow, the peak in the [O/H] distribution function is at higher 
metallicity than the [O/H] at which the galaxy is when at its peak 
star formation rate. Since M*,o ~ 10 11 M© galaxies are skirting 
the peak of the Z s -M^/Msfr relation (bottom-right panel of Fig- 
ure!]}, they have the nearly maximal stellar [a/H] distribution; the 
turnover in the metallicity scaling relations at high Kh and high 
Msfr is why the log M±/Mq = 10.5 and 11.5 distribution func- 
tions are surprisingly similar. 

Table[T]gives the distribution functions plotted in the top panel 
of Figuref2] At log M*/Mq < 9.25 (i.e., the regimes that are cal- 
culated entirely from extrapolations), the distribution functions at 
are fairly well described by power-laws of the form 



log 



AM 



Alog(0/H) 



= a[0/H] +b, 



(5) 



where AM is the mass fraction in a bin of width Alog(0/H), 
[O/H] = log[(0/H)/(0/H) ], and a and b are functions of stel- 
lar mass such that 



a = 0.2205 log M*/M Q -0.9305 and 
b = -0.42 log M*/M +4.36. 



(6) 
(7) 



At higher masses, this power-law distribution turns over, such that 
the high-metallicity distribution is well-described as a truncated 
Gaussian, where the maximum metallicity is given by the z = 
gas-phase metallicity (top-right panel of FigureQ}. 

The bottom two panels of Figure[2]show the galaxy-averaged 
a-element stellar metallicities (a*) as a function of stellar mass 
at z — (middle panel) and at 2 = 0.7 and 2 = 1.1 (bot- 
tom panel). The open symbols in the middle panel show the stel- 
lar mass-weighted mean a^-M* relation. Though it is a small ef- 
fect, mass loss must be taken into account to correctly compute 
the m ass-weighted average for a n entire galaxy; we use the meth- 
ods of lLeitner & Kravtsovl fcOllh for consistency with our adopted 
star formation histories. The solid symbols denote a B-band lu- 
minosity weight ed mean, which we obtain using the stellar popu- 
lation models of iBruzuaT & Charlotl J200I) and our assumed star 
formation history. The B-band luminosity weighted mean metal- 
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Figure 2. Top: Distribution of stellar [O/H] in galaxies of different M <: o- 
Values are given in Table[T] Measurements of the Milky Way disk from 
iBensbv et al] <2004l) are shown for reference. Dashed lines denote extrap- 
olation. Middle: Th e stellar metallicity-mass relation; data from SDSS 
iGallazzi et"ai]|2005l) is shown as the black solid line (dash ed lines show- 
ing H o~ dispersion) using the [a/Fe]-M* conversio n from [ T homas e t alj 
<2005l) . The fit to the Z+-M* relation for dwarfs bv lWooetalJ J2008t) is 
shown as a dot-dashed line. The open symbols denote mass-weighting, 
while the filled symbols are the a abundances weighted by the B-band lu- 
minosity. Bottom: The a±-M+ relation at z = (diamonds), 0.7 (circles), 
and 1.1 (triangles). Symbols of the same color denote tracks of individual 
galaxies. Values are given in Table|2] 
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licity closely appro ximates observational estim ates of metallicity 
from spectral lines dTrager & Somervilldl2009h . For comparison, 
we plot the SPSS ste llar me tallicity-mass r elation as measured by 
IGallazzi et all {2005)0 The IGallazzi et al.l measurements are pri- 
marily sensitive to iron abundances (A. Gallazzi, private communi- 
cation), but our model is implicitly tied to o xygen abundan ces (i.e., 
a elements). Therefore, we renormalize the lGallazzi et alj data us- 
ing th e [a/Fe] ratios (relative to Solar) found by Thomas et al.1 
(2005), where 



[a/Fe] = 0.16 + 0.062(log M*/M© - 10). 



(8) 



A major caveat of this renormalization is that it is based on mea- 
surements of a/Fe for elliptical galaxies, which are predominately 
not forming stars, with an intermediate step of connecting to 
galaxy velocity dispersions. Furthermore, the lGallazzi et alJd2005h 
data are for both star-forming and non-star-forming galaxies; in 
particular, the M* < 10 10 ' 5 M© steep part of the relation is pre- 
dominately star forming galaxies, while the flatter higher mass por- 
tion is predominately passive galaxies. At M*,o > 10 1 M©, the 
predictions are not robust because of extrapolation of Equation l[2} 
to very high star formation rates0 

Overall, the agreement of our simple model predictions with 
the observed average stellar metallicity relation is remarkably good 
over ~ three orders of magnitude. Our prediction for the mean 
stellar metallicity is in excellent agreement with the lGallazzi et al.1 
results at around the stellar mass of the Milky Way 
(fewxlO 10 Mq), where they are expected to be the most robust. 
At stellar masses above ~ 10 11 Mr , our p redicted metallici- 
ties turn down, while the lGallazzi et alj d2005l ) results continue to 
increase — this is likely because of the breakdown of our assump- 
tions at these high masses, where the local galaxy population is 
dominated by galaxies that are not on the star forming main se- 
quence, and which are like ly to have experienced significant growth 
throu gh gas-poor mergers dDe Lucia & Blaizotl20 07; Skelt on et al.1 
l2009h . 

Between 10 9 < M* < 10 10 M©, our predictions overshoot 
the lGallazzi et alj results slightly but are ag ain in very good agree- 
ment with the results of IWoo et alj j2008h based on local group 
galaxies from 10 8 to a fe w x 10 9 M© (a relation th at extends four 
dex further down in M*; [Kirbv et al. 2010. . 201 lab . However, it is 
importa nt to note that the | Leitnerl | |2012|) star formation histories 
(and the lMannucci et alj|20 id metallicity relation) used in our pre- 
dict ions are entirely extrapolated in this regime. At M* ~ 10 8 , 
the iLeitnerl extrapolated star formation histories differ substan- 
tially from the histories derived from the "fossil record" (resolved 
color magnitude d ata for individual stars) in nearby dwarf galaxies 
dWeisz et alj|201 ll ). The Leitner estimates would imply that dwarf 
galaxies form most of their mass at z < 1, while observat ionally, 
downsizing reverses in this regime: the lWeiszetatlbOllI) results 
imply that they form most of their stars at z > 2 (see lLeitnerll2012l 
for a detailed discussion). In our framework, however, z = metal - 
licities do not depend on when the metals were made: as long as 
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Table 2. Evolution of the stellar metallicity-mass relation, as plotted in the 
bottom panel of Figure|2] stellar masses M» are in units of M© and metal- 
licities a* are relative to »©. Rows denote the growth of individual galax- 
ies. 



the ensemble-averaged galaxies remain on the M^-MsFR-Zg rela- 
tion, the average stellar metallicities will be well-described by our 
modelEI 

Thus, at both low and high masses, it is difficult to compare 
our predictions to the current observations. Based on this combi- 
nation of observational factors and the regimes in which we are 
extrapolating to get predictions, we consider our predictions to be 
the most robust at stellar masses roughly comparable to the Milky 
Way (M*, ~ 1 to 5 x 10 10 Mq). 

Finally, the steep drop in predicted metallicities at 
logM*/M© < 9.5 is likely caused by an underestimate of 
12 + log(0/H) from extrapolating Equation© outside the 
observed range of M* and Mspr. Though the fit given by 
lLara-Lopez et al.1 (equation[4]l gives nearly the sa me a+ at low 
as shown in Figured the piecewise-linear fit from lMannucci et"al] 
(their equation 5) gives higher metallicities for low-mass galaxies 
(log a*/a© ~ -0.5 at M* ~ 10 8 M©), though it dramatically 
over-predicts the metallicities of high-mass galaxies. 

The bottom panel of Figure|2] and Table|2] give the predicted 
evolution of the a* -Mi, relation from z = 1.1, 0.7, to 0.00 The 
stellar metallicity-mass relation is slightly steeper at higher redshift 
than at z = 0; this evolution is driven by the relatively late buildup 
of stellar mass in low-mass galaxies relative to more massive galax- 
ies, as can be seen in the top panel of FigureQ] In particular, at 
Mi, > 10 10 5 M©, where ongoing surveys are measuring stellar 
metallicities at z ~ 0.7 for the first time (Gallazzi et al. in prepa- 
ration), we find a typical decrease in a* of less than ~ 0.1 dex for 
star-forming galaxies. (Note that once star formation has ceased, 
stellar metallicities will not evolve; thus if passive galaxies at z = 
have not formed stars since z — 0.7„ then galaxies will either re- 
main in place or move to the right in the plane through 
stellar mass growth via dry mergers). This rather modest evolution 



4 We note that our model predicts metallicities for stars formed in situ, 
whereas many star-forming galaxies have a bulge component that is thought 
to have been formed from the merg ing of smaller galaxies (e.g.. iToomrd 
ll977l : lKormendv & Kenn icutt 2004). Because the SDSS spectra are taken 
from fibers placed on the central 3" of each galaxy, the derived metallicities 
could potentially be highly sensitive to such a bulge component. 

5 Interestingly, there is increasing evidence that at high the depen dence 
in 12 + log(0/H) on star formation rate reverses iYates et al.l2012h . 



6 Furthermore, it is intriguing to note that the observed average stellar 
metallicities in this mass range obey such a tight relation desp ite the wide 
range in metallicity distribution functions (Dolphin et al. 2005) and the rel- 
atively wide spread in gas-phase abundances (e.g. JzahMet*al.l2012al) . 

7 Thou gh stellar metallicities h ave been measured at redshifts as high as 
z ~ 3 ISommar iva et aljl2()12l). we do not provide predictions for stellar 
metallicities at higher redshift due to uncertainties in the star formation his- 
tories. 
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is in stark contrast to the evolution seen at all stellar masses in the 
Z g -M* relation (top-right panel of Figure [T}; this difference is pri- 
marily due to the earlier buildup of stellar mass in the more massive 
galaxies. 



4 CONCLUSIONS 

We have combined empirically derived star -formation his tories 
from observations of high redshift galaxies dLeitnenl2012h with 
the empirically determined relation bet ween stellar mass, star for- 
mation rate, and gas-phase metallicity dMannucci et alj|2010l) to 
derive the distribution of stellar metallicities at z = and the 
evolution of the galaxy-averaged stellar metallicity-mass relation 
(Figure[2] and Tables[T] and [2}. We find that the hypothesis that 
the observed z = relation between stellar mass, star forma- 
tion rate, and gas-phase metallicity holds up to at least z ~ 2.5 
is consistent with stellar metallicities, within the limits of the cur- 
rent observations. This comparison, however, is hampered by both 
a lack of a well-measured M*-Msfr-Z & relation at low stellar 
masses, low star formation rates, and high star formation rates, 
and a lack of well-characterized a-element distribution functions 
in both the Milky Way and other star-forming field galaxies with 
10 s < M* < 1O 11 M0. From the modeling side, upcoming char- 
acterizations of the Mt-MsFR-Zg relation to lower stellar masses 
and across a broader range of star formation rates (B. Andrews & 
P. Martini, in preparation) will help alleviate the need for extrap- 
olation when calculating 12 + log(0/H). Upcoming surveys will 
also provide a more complete census of the stellar metallicity dis- 
tribution function for the Milky Way (e.g., the SDSS-III APOGEE 
survey) and Local Group dwarf galaxies (E. Kirby et al. in prepa- 
ration). 

We have assumed here that all of the scatter in 12 + log(0/H) 
at fixed stellar mass is related to variations in the star formation 
rate, but the resulting relation still has scatter that is cor related 
with, e.g., galaxy size tellison et alj2008al:l Yabe et al.l2012l) or en- 
vironment jKewlev et alj2006l : Ellison et al.ll2008bl ; |Pasquali et al .1 
120121 : lHughesetal.ll2012h . Moreover, the scatter i n Z+ at fixed 
stella r mass is known to correlate with galaxy age dGallazzi et al .1 
120051) . In our framework, this can be seen as variations in galax- 
ies' star formation histories via t he ~ 0.3 dex scatter in star for- 
mation rate at fixed stellar mass l lNoeskeetal1l2007l ;l Karim et al .1 
l201lHLeitnenl2012l) . Thi s variation can be p artially attributed to 
changes in environment jPasquali et alj|201fj|) . especially if envi- 
ronm ent introduces additi onal scatter to the Z g ~M*-MsFR relation 
(e.g JPasquali et alj|2012l) . 

We have also assumed that the buildup of all of the stellar mass 
in a galaxy can be well described by Equation (TJl. This is clearly 
not the case, as, e.g., there are stars in the Milky Way that were 
formed before z ~ 2.5 and with log Z*/Zq < —3. However, these 
stars comprise a very small fraction of the total stellar mass or light, 
at least in the Milky Way, and could be addressed by more complete 
model of star formation histories, including minor mergers which 
bring in a populations of lower-metallicity stars. 

One of the strengths of the method presented here is that 
it does not require us to explain why the relationship between 
gas-phase metallicity, stellar mass, and star formation rate appar- 
ently does not evolve with redshift by invoking carefully evolv- 
ing gas fractions, accretion rates, o r outflow efficien c ies (such as 
in, e. g., IPeeples & Shankail 1201 ll ; iDave et al] |2012| ; iDaval etHll 
120121) . On the other hand, integrating over the evolution of the 
mass-metallicity relation does neatly predict what fraction of Oxy- 
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gen ever produced by a z — galaxy should be still locked up 
in stars, and how steeply this fraction should increa s e with stel- 
lar m ass (c.f.. lGallazzi et alj|2008l ; iKirbv et al.ll201 lbl : IZahid et al.l 
120121) . Moreover, by combining this kind of calculation with mea- 
surements of the mass-metallicity r elation and z = gas frac- 
tions (e.g.. |Peeples & Shankar 201lh. one can place constraints on 
the to tal amount of Oxygen galaxies have expelled over their life - 
time jBouche et all 12001 IKirbv et alj |2011bl : IZahid et al.l |2012|) . 
though for star-forming galaxies, uncertainties in the gas-phase 
abundances and relevant gas masses make this calculation difficult. 
More globally, it is especially intriguing that the implied masses of 
both metals and all baryons expelled by galaxies over their lifetime 
appears potentially consistent with the oxygen mass observed in 
the ci rcumgalactic medium of z ~ 0.25 galaxies dTumlinson et al.l 
I20T l|; J. Werk et al. in preparation). 
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